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Abstract 

Environmental factors such as diet, intake of vitamin D supplements and exposure to sunliglit are l<nown to influence serum 
vitamin D concentrations. Genetic epidemiology of vitamin D is in its infancy and a better understanding on how genetic 
variation influences vitamin D concentration is needed. We aimed to analyse previously reported vitamin D-related 
polymorphisms in relation to serum 25{OH)D concentrations in 201 healthy Danish families with dependent children in late 
summer in Denmark. Serum 25(OH)D concentrations and a total of 25 SNPs in GC, VDR, CYPIRI, CYP24A1, CYP27B1, C10or88 
and DHCR7/NADSYN1 genes were analysed in 758 participants. Genotype distributions were in Hardy-Weinberg equilibrium 
for the adult population for all the studied polymorphisms. Four SNPs in CYP2R1 (rsl 562902, rs71 16978, rsl 0741 657 and 
rsl0766197) and six SNPs in GC (rs4588, rs842999, rs2282679, rsl2512631, rsl6846876 and rsl7467825) were statistically 
significantly associated with serum 25{OH)D concentrations in children, adults and all combined. Several of the SNPs were in 
strong linkage disequilibrium, and the associations were driven by CyP2/?7-rsl0741657 and rsl0766197, and by GC-rs4588 
and rs842999. Genetic risk score analysis showed that carriers with no risk alleles of CyP2R7-rsl 0741 657 and rsl 07661 97, 
and/or GC rs4588 and rs842999 had significantly higher serum 25(OH)D concentrations compared to carriers of all risk 
alleles. To conclude, our results provide supporting evidence that common polymorphisms in GC and CYP2R1 are associated 
with serum 25(OH)D concentrations in the Caucasian population and that certain haplotypes may predispose to lower 
25(OH)D concentrations in late summer in Denmark. 
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Introduction 

Vitamin D deficiency is a widespread problem in developed 
countries [1]. Severe vitamin D deficiency causes osteomalacia, or 
childhood rickets, osteoporosis and fractures because of reduced 
calcium absorption [2] . Low vitamin D concentrations may also be 
related to various non-skeletal health outcomes, including cardio- 
vascular diseases [3] , obesity [4] , diabetes [5] , asthma [6] , multiple 
sclerosis [7] , occurrence of a large range of cancer diseases [8] and 
overall mortality [9,10]. 

In humans, vitamin D is produced mainly in the skin during 
exposure to solar ultraviolet blue (UVB) radiation (270-300 nm) 
[11]. UVB radiation converts 7-dehydrocholesterol (7-DHC) in 
the skin to pre-vitamin D3, which immediately undergoes a 
thermal isomerization to vitamin D3. Dietary sources provide two 
forms of vitamin D: Vitamin D2 (ergocalciferol) derived from 
invertebrates (plants and fungi) and vitamin D3 (cholecalciferol) 
derived from animal sources. Ingested vitamins Dj and D3 are 
absorbed in the small intestine and transported with chylomicrons 
and lipoproteins to the liver, whereas dermally synthesized vitamin 



Di; diffuses via the blood to the liver tightly bound to group- 
specific complement (GC) [12]. 

Dietary or dermally synthesized vitamin D (hereafter "D" refers 
to D2 and D3) undergoes a series of enzymatic conversions in the 
liver and kidneys to become biologically active. The hepatic 
enzyme 25-hydroxylase (CYP2RI) converts vitamin D to 25- 
hydroxyvitamin D (25(OH)D). This is the major circulating form 
of vitamin D in the blood. To become biologically active, 
25(OH)D is converted to 1,25-dihydroxyvitamin D 
(l,25(OH)2D). This occurs mainly in the kidneys, but also in 
other tissues expressing the enzyme 25(OH)D-lot-hydroxylase 
(CYP27B1). The biological effect of vitamin D is mediated when 
l,25(OH)2D binds to the vitamin D receptor (VDR). To prevent 
excessive vitamin D signalling in the target organs, l,25(OH)2D 
limits its own activity by inducing 24-hydroxylase (CYP24A1) 
converting l,25(OH)2D to the biologically inactive water-soluble 
calcitroic acid which is excreted in the bile [1,12,13]. 

The best biomarker of vitamin D concentration is the serum 
25(OH)D concentration. Approximately 25% of the inter-individ- 
ual variability in plasma 25(OH)D concentrations can be 
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explained by external factors such as diet, regular use of vitamin D 
supplements and exposure to sunlight (dependent on season and 
latitude) [14,15]. Genetic factors may contribute to vitamin D 
concentrations. Results from twin and family-based studies 
indicate that blood vitamin D concentrations to some extent are 
under genetic control. The results have been inconsistent with a 
wide variabihty in heritabUity estimates ranging from 23 to 80% 
[15-21]. Furthermore, ethnic differences in vitamin D concentra- 
tions have also been described [22]. 

Genetic epidemiology of vitamin D is in its infancy and a better 
understanding of how genetic variation influences vitamin D 
concentrations is needed. A growing number of studies have 
uncovered polymorphisms associated with vitamin D concentra- 
tions. By candidate gene analysis, five genes have been found, 
including GC, CrP24Al, CrP2Rl, CrP27Bl and VDR [23]. 
Recently, two genome- wide association studies (GWAS) of vitamin 
D [24,25] confirmed the associations of common variants in GC 
and CYP2R1 genes. Furthermore, nicotinamide adenine dinucle- 
otide synthetase- 1/7-dehydrocholesterol reductase (NADSYNl/ 
DHCR7), and the region harbouring the open-reading frame 88 
(C10orf88) on chromosome 10q26.I3 were also found to be 
associated with vitamin D concentrations in blood. 

In Denmark, low vitamin D status is common during the winter 
due to inadequate dietary intakes and lack of solar radiation from 
September to April [26]. We assessed vitamin D status in late 
summer (September to October), where the Danes vitamin D 
concentration peaks but are not saturated [27], in families with a 
broad span in age in both children and adults. In children, the role 
of genetic variation in determining serum 25(OH)D concentra- 
tions is an understudied area. 

In this study, we analysed previously reported vitamin D-related 
polymorphisms in relation to serum 25(OH)D concentrations in 
201 healthy Danish families with dependent children to confirm 
previous findings and thus help identifying individuals that may 
have increased risk of developing vitamin D insufiSciency. 

Subjects and Methods 

Study population 

The present cross-sectioned study used baseline data from the 
VitmaD intervention study described in detail elsewhere [28]. 
Briefly, 201 Danish families with dependent children (n = 782) 
were enrolled. The partidpants w("re 4- to 6()-years old. Baseline 
blood samples were collected in September and October 20 1 0 and 
were obtained from 770 participants. The study was conducted 
according to the guidelines in the Declaration of Helsinki and the 
protocol was approved by the Research Ethics Committee of the 
Capital Region of Denmark (H-4-20 10-020) and registered at 
http://clinicaltrials.gov (NCTOl 184716). All adult participants 
and guardians on the behalf of the children participants gave 
written consent to participate. 

DNA extraction and genotyping 

DNA was extracted from peripheral blood leukocytes as 
described by Miller a al. [29] and stored in TE-bufier at -80°C. 
The DNA was diluted to 10 ng/jil using a Nanodrop® ND-1000 
Spectrophotometer (Thermo Fisher Scientific Inc., Wilmington). 
Single nucleotide polymorphisms (SNPs) were genotyped using the 
Sequenom MassARRAY iPLEX Gold platform (Sequenom, San 
Diego, California) at the Department of Biomedicine, Aarhus 
Universit\r, Denmark. Genotyping was successful for 762 partic- 
ipants (99.0'/o). To confirm the accuracy of genotyping duplicate 
samples (10%) yielded 100% reproducibility. 



All SNPs were located in or near genes involved in vitamin D 
synthesis, activation or degradation. The following SNPs were 
selected on the basis of evidence of significant association in 
previous studies: CYP2R1 (rsl562902; rs71 16978; rsl0741657; 
rsl0766197) CYP24A1 (rs229624; rs2426496; rs4809960; 
rs60 13897; rs 1 7219315) CYP27B1 (rs 108770 12) ClOorfSS 
(rs6599638) DHCR7/NADSYN1 (rsl790349; rsl2785878) GC 
(rs4588; rs222020; rs842999-triaUeUc; rs2882679; rs2298849; 
rsl2512631; rsl6846876; rsl7467825) VDR (rs731236 (TaqI), 
rs757343 (Trul); rs7139166; rsl0783219). 

Deviation from Hardy-Weinberg equilibrium (HWE) was 
tested for the adult population using Chi-square test with 
Bonferroni's correction (P-value 0.05/25 SNPs = 0.002). No 
significant deviation from HWE was observed. Linkage disequi- 
librium (LD) between polymorphisms was evaluated using 
Pearsons' r, SNAP version 2.2 (http://www.broadinstitute.org/ 
mpg/snap/ldsearchpw.php) and Haploview software version 4.2 
for the adult population. 

Measurement of serunn 25(OH)D concentrations 

Measurements of serum 25(OH)D concentrations are described 
in detail elsewhere [28]. Briefly, blood samples were obtained 
without prior fasting and serum was stored in aliquots at -80°C 
until analysis. Measurements of serum 25(OH)D concentrations 
relied on the determination of both 25(OH)D2 and 25(OH)D3 and 
were conducted by isotope dilution liquid chromatography 
tandem mass spectrometry (LC-MS/MS) at Clinical Biochemical 
Department, Holbiek Hospital, Denmark. As primary calibrator 
the standard reference material, vitamin D in humans (SRM 972) 
from the National Institute of Standards and Technology was 
used. The analytic quality of 25(OH)D assay was assured by 
Vitamin D External Quality Assessment Scheme certification and 
the mean bias was -3.2%. The Inter-assay CVs for 25(OH)D2 
were 7.6'/o and 4.6% at 43 and 150 nmol/L, respectively, and for 
25(OH)D:5 2.2% and 2.8% at 30 and 180 nmol/L, respectively 
[28]. Of the 762 participants that were successfully genotyped, 
baseline serum 25(OH)D concentrations were measured for 758 
participants. 

Statistical analysis 

Statistical analyses were performed using SAS Enterprise Guide 
4.3 (SAS Institute, Inc., Gary. USA). Serum 25(OH)D concentra- 
tions were log transformed to approximate a normal distribution 
and all means are presented as geometric means. A nominal P- 
value of 0.05 was considered statistically significant. Linear mixed 
models with family as a random factor were applied to account for 
the possible dependence b(;tween the participants. Furthermore, in 
the linear mixed models the following categorical variables were 
included: age (4-11, 12-17, 18-40, 41-60 years), sex (male, 
female), BMI (underweight, normal weight, overweight, obese) 
according to standards for children [30] and the WHO 
International standards for adults [31], ski or sun holidays (yes, 
no), solarium use at least once a week (yes, no), dietary vitamin D 
(quartiles: <1.7, 1.7-2.4, 2.5-3.3 and >3.3 Hg/d), multivitamin 
and vitamin D supplement users (yes, no). The data were obtained 
from a self-administered web-based questionnaire and a semi- 
quantitative food frcfjuc-nc)' rjuestionnaire based on the last six 
months. Pearson's r were calculated on the adult population and 
were used to assess the degree of hnkage between Unked SNPs. 
Haplotypes were inferred manually among the adults, only since 
the children were not population-based. The inferred haplotype 
combinations described lOO'/o and 97% of the observed genotypes 
among the adults for CTP2R1 and GC genes, respectively. Among 
the children the inferred haplotype combinations described 100% 
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and 96% of the observed genotypes for CYP2R1 and GC genes, 
respectively. Each derived haplotype was assigned a number. 
Homozygote haplotype combinations were numbered with two 
identical numbers e.g. 11. The combinations of heterozygote 
haplotypes were given by the combination of the number of each 
haplotype e.g. 1+ 2 = 12. 

Genetic risk scores were calculated as the sum of risk alleles and 
included as risk factors in linear mixed models adjusted for family 
and confounding variables. The correlation coefficient for 
rsl0741657, rsl0766197, rs4588 and rs842999 were very similar 
and therefore it was not necessary to weight the score by elfect 
size. AH the analyses were performed separately for children, 
adults and for all combined. 

Results 

Genotyping and serum 25(OH)D concentrations were available 
for 758 participants. Table 1 summarizes the basic characteristics 
of the study population, previously described in detail elsewhere 
[28]. The median age among children was 10 years (range: 4 to 17) 
among adults 41 years (range: 18 to 60) and for all combined 30 
years. 

Associations between genotypes and serum 25(OH)D concen- 
trations are shown for children, adults and aU combined in Table 
2. After adjustment for family and confounding factors, all four 
analysed SNPs in CTP2R1 were statistically significantly associated 
with serum 25(OH)D concentrations in all three groups. Further- 
more, for all three groups none of the analysed SNPs in CTP24A1, 
CrP27Bl, ClOorJSS and DHCR7 /NADSYNl were statisticaUy 
significantly associated with serum 25(OH)D concentration. For 
all three groups all analysed SNPs in GC, except rs2298849 (in all 
three groups) and rs222020 (in adults and all), were statistically 
significantiy associated with serum 25(OH)D concentration. The 
VDR rs731236 was only statistically significantly associated with 
25(OH)D concentration in all combined and rs757343 was 
statistically significant in children and all combined. Only SNPs 
that were statistically significantly associated with 25(OH)D 
concentrations in children, adults and all combined were included 
in further analyses. 

Haplotype and genetic risk score analysis of CYP2R1 

In the adult population, rsl0741657-rs7 1 16978 (Pearson's r 
= 0.90), and rsl076697-rsl562902 (Pearson's r =-0.86, data not 
shown) were in strong LD. To establish which of the SNPs had the 
strongest association to serum 25(OH)D concentrations, we assess 



the association between one SNP and serum 25(OH)D concen- 
trations while adjusting for the other SNPs, family and confound- 
ing factors in a linear mixed model. After adjustment, rsl0766197 
(p = 0.0846) had the strongest association compared to rs 1562902 
(p = 0.8211), and rsl0741657 (p = 0.2545) had the strongest 
association compared to rs7 116978 (p = 0.3087, data not shown). 
In further analysis only rsl076697 and rsl0741657 were included. 

The two CrP2Rl variants rsl0741657 and rs7 116978 formed 
four haplotypes, where haplotype 1 and 2 were most frequent 
(Table 3). The possible combinations of the four homozygote 
haplotype are shown in table 3. One genotype combination could 
be assigned to both haplotype combinations 1 2 or 34, but based on 
the observed haplotype frequencies, the most likely combination 
was 12. After adjustment for family and confounding factors, 
carriers of 2 copies of the AG-haplotype (haplotype combination 
33) had the highest mean serum 25(OH)D concentration (73.8 
(60.1-90.6), 72.9 (57.3-92.5) and 81.3 (66.4-99.6) nmol/L) in 
children, adults and all combined, respectively. In a linear mixed 
model, only the homozygous haplotype combinations were 
included and haplotype combination 44 was excluded because 
only two participants carried this haplotype combination. The 
homozygous haplotype combinations were significandy associated 
with serum 25(OH)D concentrations (p = 0.0059, 0.0450 and 
0.0007) in children, adults and all combined, respectively. 

We calculated a genetic risk score (range 0-4) as the sum of the 
number of G-alleles of rsl0741657 and A-aUeles of rsl0766197 
(Figure 1, A). After adjustment for family and confounding factors, 
carriers of no risk alleles had significantly higher serum 25(OH)D 
concentrations (74.0 (60.3-90.0), 73.0 (57.5-92.6) and 81.3 (66.4- 
99.5) nmol/L) compared to carriers of all four risk alleles (61.2 
(57.5-92.6), 64.0 (50.6-80.9) and 69.8 (57.0-85.4) nmol/L) in 
children, adults and all combined, respectively. Overall, there was 
20.9, 14.1 and 16.5% difference in serum 25(OH)D concentra- 
tions between carrying no risk alleles and carrying all four risk 
alleles in children, adults or all combined, respectively. 

Haplotype and genetic risk score analysis of GC 

In the adult population, rs4588 was in strong LD with 
rs2282679 (Pearson's r = 0.997), rsl7467825 (Pearson's r = 0.997) 
and rsl6846876 (Pearson's r = 0.805). Furthermore, rsl7467825- 
rs2282679 (Pearson's r= 1.00), and rs2282679-rsl6846876 (Pear- 
son's r =0.8021, data not shown) were also in strong LD. To 
establish which of the 4 SNPs had the strongest association to 
serum 25(OH)D concentrations, we assess the association between 
one SNP and serum 25(OH)D concentrations while adjusting for 



Table 1. Basic characteristics of the study population and determinants of serunn 25(OH)D concentrations 





Characteristics 


Children 


Adults 


All 


Number 


348 


414 


762 


Female/Male (n/n) 


181/167 


209/205 


390/372 


Age, median (range) 


10 (4-17) 


41 (18-60) 


30 (4-60) 


BMI (kg/m^)» 


17.44±2.89 


25.47±4.30 


21.79±5.45 


Serum 25(OH)D (nmol/L)* 


74.38±17.31 


74.87±21.70 


74.65 ±19.82 


Dietary Vitamin D (jjg/d)* 


2.69±1.35 


2.96 ±2.04 


2.84±1.77 


Multivitamin or vitamin D supplement users (yes/no) 


141/203 


1 1 3/297 


254/500 


Solarium use (yes/no) 


2/342 


10/401 


12/743 


ski or sun holidays (yes/no) 


195/149 


220/191 


415/340 



•Mean ± SD. 

doi:l 0.1 371 /journal.pone.0089907.t001 
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a. CYP2R1 (rsl0741657 and rsl0766197) 



b. SC(rs4588andrs842999) 
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c. CYP2R1 (rsl0741657 and rsl0766197) and GC (rs4588 and rs842999) 
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Figure 1. Genetic risl( score for CyP2R7 (rs10741657 and rs10766197) (Figure A), Cr(r45S8 and rs842999) (figure B) and CyP2R7 
(rs1 0741 657 and rs1 07661 97) and CC(r4588 and rs842999) (figure C) in cKiiidren, adults and all combined. X-axis stands for the sum of 

risk alleles. Y-axis stand for serum 25(OH)D (nmol/L). Errors bars stand for 95%-confidence interval and serum 25(OH)D concentrations are given as 
geometric means. Linear mixed models with family as a random factor, adjusted for age, sex, BMI, ski and sun holidays, solarium use at least once a 
week, dietary vitamin D intake, multivitamin and vitamin D supplement users was conducted to compare sum of risk alleles and serum 25(OH)D 
concentrations. Increasing number of risk alleles give rise to decreasing 25(OH)D concentrations. 
doi:1 0.1 371 /journal.pone.0089907.g001 



the other SNPs, family and confounding factors in a hnear mixed 
model The strongest association was observed for rs4588 
(p = 0.0099) compared to rs2282679 (p = 0.0230), rsl 7467825 
(p = 0.0230) and rsl 6846876 (p = 0.5669, data not shown). Further 
analyses only included rs4588. None of the other GC-variants were 
in LD. 

The three significant GC-variants (rs4588, rs842999, and 
rsl2512631) formed five haplotypes, where haplotype 1 and 2 
were the most frequent (Table 4). The combinations of the five 
haplotypes are shown in table 4. The five haplotypes could explain 
723 of the 762 (95%) observed genotype combinations in GC (data 
not shown). The association between haplotype combinations and 
serum 25(OH)D concentrations was statistically significant in 
children (p = 0.0344), and all combined (p = 0.0018) but not in 
adults (p = 0.1541). 

Carriers of haplotype combination 22 encompassing the variant 
alleles of rs4588 and rs842999 had low serum 25(OH)D 
concentrations. Conversely, carriers of haplotype combination 
1 1 encompassing the variant allele of rsl 25 12631 had high serum 
25(OH)D concentration. Thus, the variant allele of rsl 25 12631 
was associated with high low serum 25(OH)D concentrations and 
the variant alleles of rs4588 and rs842999 were associated with low 
serum 25(OH)D concentrations. Since the lowest serum 25(OH)D 
concentrations were observed for haplotype combination 22 



carriers, this could indicate that rs4588 is the biologically relevant 
polymorphism rather than rs842999 since haplotype combination 
44 encompassing the C-aUele of rs842999 is associated with higher 
serum 25(OH)D concentrations. 

The genetic risk score (range 0-4) was calculated as the sum of 
the number of A-alleles of rs4588 and C/A-alleles of rs842999 
(Figure 1, B). After adjustment for famUy and confounding factors, 
we found that an increasing number of risk alleles was associated 
with lower serum 25(OH)D concentrations. Carriers of no risk 
alleles had significantly higher serum 25(OH)D concentrations 
(68.1 (56.2-82.6), 81.0 (64.2-102.2) and 86.5 (70.9-105.5) nmol/ 
L) compared to carriers of all four risk alleles (50.3 (40.3-62.7), 
67.5 (53.6-84.9) and 70.1 (57.2-84.8) nmol/L) in both children, 
adults and all combined, respectively. Overall, there was a mean 
difierence in 25(OH)D concentrations of 35.4, 20.0 and 23.4% 
between carrying no risk alleles and carrying all four risk alleles in 
children, adults and all combined, respectively. 

For the tri-aUelic variant rs842999, there was a dose-dependent 
relationship between serum 25(OH)D concentrations and carriage 
of none, one or two copies of the G-allele (Figure 2). Thus, carriers 
of two copies of the G-allele, had statistically significantly higher 
serum 25(OH)D concentrations (69.2 (56.8-84.3), 79.0 (62.8-99.4) 
and 84.8 (69.6-103.4) nmol/L) compared to carriers of only one 
G-allele (65.6 (53.9-79.9), 73.7 (58.8-92.4) and 79.0 (64.9-96.1) 
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Figure 2. Dose-dependent relationship between genotype GG, GX and XX of rs842999 and serum 25(OH)D concentrations. X-axis 
stands for genotype GG (GG), GX (GC or GA) and XX (CC, CA or AA) of rs842999. Y-axis stand for serum 25(OH)D (nmol/L). Errors bars stand for 95%- 
confidence interval and serum 25(OH)D concentrations are given as geometric means. Linear mixed models with family as a random factor, adjusted 
for age, sex, BMI, ski and sun holidays, solarium use at least once a week, dietary vitamin D intake, multivitamin and vitamin D supplement users was 
conducted to compare rs842999 genotypes with serum 25(OH)D concentrations. There was a dose-dependent relationship between serum 25(OH)D 
concentrations and carriers of none, one or two copies of the G-allele. Carriers of two copies of the G-allele, had higher serum 25(OH)D 
concentrations compared to carriers with only one G-allele or non-carriers in children, adults and all combined, respectively. 
doi:10.1371/journal.pone.0089907.g002 



nmol/L) in children, adults and all combined, respectively. The 
lowest serum 25(OH)D concentrations were observed in non- 
carriers of the G-aUele (59.5 (48.7-72.6), 67.4 (53.8-84.4) and 72.8 
(59.7-88.8) nmol/L) in both children, adults and all combined, 
respectively. 

Finally, we made a joint genetic risk score analysis including 
CrP2Rl (rsl0741657 and rsl0766197) and GC (rs4588 and 
rs842999) (Figwe 1, C). The genetic risk score (range 0-8) was 
calculated as the sum of the number of G-alleles of rs 1074 165 7, A- 
alleles of rsl0766197, A-alleles of rs4588 and C/A-alleles of 
rs842999 (Figure 1, C). After adjustment for family and 
confounding factors, carriers of no risk alleles had statistically 
significantiy higher 25(OH)D concentrations (78.4 (63.6-96.7), 
86.3 (66.1-112.7) and 89.0 (72.0-110.0) nmol/L) compared to 
carriers of all eight risk alleles 43.4 (32.4-58.2), 55.3 (37.5-81.4) 
and 53.0 (39.6-70.9) nmol/L) in children, adults and aU 
combined, respectively. Overall there was a mean difference in 
25(OH)D concentrations of 80.6, 56.1 and 67.9% between 
carriage of no risk alleles and carriage of all four risk alleles in 
children, adults and all combined, respectively. 

Discussion 

In this present study, we studied the association of 7 prominent 
vitamin D-related genes with serum 25(OH)D concentrations in 
201 Danish families with dependent children in late summer in 
Denmark, and found that common variants in CTP2R1 and GC 
genes were statistically significantly associated with serum 
25(OH)D concentrations. 

The CTP2R1 gene encodes the key enzyme that converts 
vitamin D to 25(OH)D in the liver [12] and thus genetic variation 
in this gene might affect 25(OH)D synthesis. We found that 
CrP2Rl variants rsl562902, rs71 16978, rsl0741657 and 
rsl0766197, were significantly associated with serum 25(OH)D 
concentrations in both children, adults and all combined. 
Furthermore, rsl0741657-rs71 16978, and rsl0766197-rsl562902 
were in strong LD. The association appeared to be driven by 
rsl0741657 and rsl0766197, which are located in the promoter 
region of the CYP2R1 gene. We found that non-carriers of 



rsl0741657 and rsl0766197 risk alleles had the highest mean 
serum 25(OH)D concentrations. 

Our results are consistent with previous findings. In the study of 
Wjst et al. [21], rsl0766197 was significantiy associated with 
25(OH)D concentrations in 872 subjects from the German 
Asthma Family Study. Ramos-Lopez et al. [32] found a 
statistically significant association between rsl0741657 and serum 
25(OH)D concentrations in 203 German diabetes families. Two 
genome-wide association studies (GWAS) of vitamin D concen- 
trations were published in 2010 [24,25]. Ahn et al. [24] performed 
a combined meta-analyses in 4,501 subjects from five adult 
Caucasian cohorts and found that rs2060793, which is in LD with 
rsl0741657 (D= 1, r^= 1, HapMap Data Rel 24/phase II Nov 
08), was associated with serum 25(OH)D concentrations. Further- 
more, these findings were successfully replicated in 2,221 subjects. 
Wang et al. [25] found that rsl0741657 was significantiy 
associated with 25(OH)D concentrations in 30,000 subjects of 
European descent from 15 cohorts. In the study of Bu et al. [33], 
rsl0741657 and rsl0766197 were found to be significantiy 
associated with serum 25(OH)D concentrations in 496 unrelated 
healthy Caucasian subjects. Lasky-Su et al. [34] conducted a 
combined analysis in 1,164 subjects from two cohorts of Caucasian 
and Costa Rica asthmatic children and found that rsl0741657 was 
significantiy associated with 25(OH)D concentrations. Zhang et al, 
[35] found that rsl0766197 was significantiy associated with 
25(OH)D concentrations in 2,897 unrelated healthy Chinese 
subjects from the Shanghai Osteoporosis Study. In the study of 
Engelman et al. [36], rs2060793 (in LD with rsl0741657 as 
mentioned previously) was significantly associated with 25(OH)D 
concentrations in 1,204 women of European descent from the 
Women's Health Initiative Observational Study. AU the afore- 
mentioned studies demonstrate that variants in the CYP2R1 gene 
predicts 25(OH)D concentrations. 

The GC gene encodes the vitamin D binding protein (DBF) that 
binds and transports blood 25(OH)D and other vitamin D 
metabolites to their target organs. Less than 0.04% of blood 
25(OH)D circulates in free form (bioavaUable). Most is bound with 
high affinity to DBF (83-85%) and with lower affinity to albumin 
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(12-15%) [37]. Variants in the GC gene may afTect the DBF 
binding and bioavailability of 25(OH)D and other vitamin D 
metabolites. Thus, there may be a relationship between phenotype 
and blood 25(OH)D concentrations. 

There is accumulating evidence that variants in the GC gene are 
associated with 25(OH)D concentrations. The most studied GC- 
variants are rs4588 and rs7041, giving three common GC- 
isoforms, GCIF (rs7041-T, rs4588-C), GClS (rs7041-G, rs4588- 
C), and GC2 (rs7041-T, rs4588-A), which dilTer by amino acid 
substitutions and/or by giycosylation (Gozdzik et al. 201 1). Several 
studies have shown that vitamin D status difiers significantly 
depending on rs4588 and/or rs7041 genotype, where the A-allele 
of rs4588 and the T-allele of rs7041 are consistendy associated 
with lower 25(OH)D concentrations [17,38-45]. In agreement, we 
found that the A-allele of rs4588 is associated with lower 25(OH)D 
concentrations. There is biological support that the affmity of both 
25(OH)D and l,25(OH)2D is higher for the C-allele of rs4588 than 
for the A-allele [46]. Based on giycosylation patterns, it is 
suggested that GC2 phenotypes that is associated with low vitamin 
D concentrations should be metabolized faster. Kawakami et al. 
observed that the metabolic rate was indeed higher in GC2-2 
individuals than in GCl-1 individuals [47]. In addition, the GC2 
genotype, which is associated with low 25(OH)D concentrations, is 
also associated with low mean DBF [43]. Strangely, the GC2 
genotype is more frequent in populations living in northern 
climates [48]. 

Since the two GWAS studies [24,25] found a strong association 
between rs2282679 and 25(OH)D concentrations, there has been 
increased focus on this polymorphism. Several studies have been 
pubUshed supporting the finding [22,34,35,49-51]. The GWAS 
GC variant rs2282679 is in high LD widi rs4588. Wang et al. [25] 
did not include rs4588 because it is not in the HapMap dataset. In 
one study sample the authors found that rs4588 was in LD with 
several associated variants from the GWAS study. In the study of 
Lu et al. [45], rs4588 and rs2282679 (r^ = 0.97) were significandy 
associated with 25(OH)D concentrations in 3,210 Han Chinese. In 
the study by Berry et al. [52], rs4588 was in strong LD with 
rs228697 (r^ = 0.98), and rs4588 was significandy associated with 
25(OH)D concentrations in 6,551 subjects from the British birth 
cohort. Zhang et al. [35] found that 2282679 and rs4588 were in 
strong LD in 2,897 unrelated healthy Chinese subjects and the 
strongest association was observed for rs4588, which accounted for 
0.7% of the variation in serum 25(OH)D concentrations. Our 
results support that rs228697 is in strong LD with rs4588 
(Pearson's r = 0.997, SNAF proxy D' = 1 r^ = 0.98) and that the 
association with serum 25(OH)D concentrations is most likely 
driven by rs4588. Zhang et al. [35] argued that it is unlikely that 
rs2282679 in itself is the disease-causing variant. The possible 
causal variant is the non-synonymous rs4588, where the C/A base 
pair change in codon 436 (previously known as 420 [36]) causes a 
Thr to Lys amino acid substitution. In agreement with Zhang et al. 
[35] we found that rs4588 was tiic strongest independent predictor 
of 25(OH)D concentrations compared to rs2282679. Furthermore, 
Zang et al [35] found that both the minor T-allele of rs4588 and 
G- allele of rs2282679 were associated with reduced DBF 
concentrations. Farticipants with 3 or 4 risk alleles of the two 
variants were more likely to have vitamin D concentrations lower 
than 50 nmol/L (20 ng/ mL) compared with non-carriers of the 
risk alleles. 

In our study, several of the significant GC variants were in strong 
LD and the strongest associations with serum 25(OH)D concen- 
trations were obser\'ed for rs4588 and rs842999. We obserx^ed a 
dose-dependent relationship between carrying none, one or two 
copies of the G-allele of the tri-alleHc rs842999 and 25(OH)D 



concentrations. Furthermore, genetic risk score analysis for rs4588 
and rs842999 showed that non-carriers of the risk alleles of rs4588 
and rs842999 had the highest serum 25(OH)D concentrations. 

We made a joint genetic risk score analysis for all four risk 
variants (CrP2i?;-rsl0741657 and rsl0766197, and GC-rs4588 
and rs842999), and found the largest% -range in mean serum 
25(OH)D concentrations (80.6, 56.1 and 67.9%) compared to 
genetic risk score analysis of CTP2R1 (rs 1074 1657 and 
rsl0766197; 20.9, 14.1 and 16.5%) or GC (rs4588 and rs842999; 
35.4, 20.0 and 23.4'M)) indicating an additive effect. In general, 
there was a better association between genetic risk score and serum 
25(OH)D concentrations in children than in adults. We speculate 
that the more risk alleles in CTP2R1 and GC genes a subject 
carries, the more prone the subject will be for having a low serum 
25(OH)D concentration. In Denmark, sufficient serum 25(OH)D 
concentrations are defined as >50 nmol/L [53]. Notably, in late 
summer in Denmark, where vitamin D status peaks in Danes, 
children carrying 7 or 8 risk alleles had insufficient serum 
25(OH)D concentrations (49.4 and 43.4 nmol/L). 

In our study population, none of the investigated SNFs in 
CrP24Al, CrP27Bl, C10orf88 or DHCR7 / NADSYNl were associ- 
ated with serum 25(OH)D concentrations. Furthermore, VDR- 
rs731236 was only statistically significant in all combined and 
rs757343 was statistically significant in children and all combined. 
False-positive (type 1 errors) results, which are common in studies 
of the association between genetic markers and outcomes, and the 
relative small sample size, resulting in statistical reduced power 
might explain these findings. We consider children and adults as 
two natural subpopulations due to biological differences, difference 
in lifestyle, eating patterns and use of multivitamins [28] . We did 
not use Bonferroni-corrected F-values because a statistically 
significant association both in children and in adults by itself 
may be considered a confirmation of an association. A limitation 
of the study is that the participants' general vitamin D status relies 
on a single measurement of serum 25(OH)D concentration. We 
were not able to calculate the genetic contribution due to the 
familiar design used in the linear mixed model. A strength of this 
study is that it is conducted in a healthy Caucasian population and 
thus the potential impact of diseases is minimized. Furthermore, 
the blood samples were collected in a relatively small geographical 
area in Denmark in September to October 2010 and analysed in a 
single batch with LC-MS/MS with low variation. Furthermore, 
many known predictors of serum 25(OH)D concentrations were 
assessed by questionnaire data. 

Genetic variants may accelerate or protect against vitamin D 
deficiency and the genetic effect is life-long. We speculate that 
individuals with genetii:ally determined low vitamin D concentra- 
tions may need different health recommendations in order to 
improve their serum 25(OI-I)D concentrations thereby avoiding 
adverse health outcomes. A study by Engelman et al. [36] found 
that in women with no risk alleles of rs4588 and rs2060793 (in 
strong LD with rs 1074 165 7 as mentioned previously) who 
consumed at least 670 lU/d vitamin D all (100%) had 25(OH)D 
> 50 nmol/L. For women carrying 1, 2 or 3-4 risk alleles and 
consuming at least 670 lU/ d vitamin D, only 84, 72, and 62% had 
25(OH)D > 50 nmol/L. Furthermore, the percentage of women 
with adequate 25(OH)D concentrations rose with each increasing 
quartile of vitamin D intake. Thus, subjects with genetic 
predisposition scc-m to benefit from dietary vitamin D supplemen- 
tation. In the study by Madsen et al. [28], vitamin Da-fortification 
of bread and mflk reduced the decrease in serum 25(OH)D 
concentrations seen during winter and ensured 25(OH)D 
>50 nmol/L in healthy Danish families. Whether such a dietary 
intervention program could ensure adequate serum 25(OH)D 
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concentrations in subjects widi genetic predisposition for vitamin 
D deficiency warrants further study. 

Conclusions 

In conclusion, our results support the current evidence that 
common genetic variation in GC and CTP2R1 may contribute to 
the variation of serum 25(OH)D concentrations in a healthy 
population. Notably, genetic risk score analysis revealed that non- 
carriers of risk alleles of CYP2R1 rsl0741657 and rsl0766197, 
and/or GC rs4588 and rs842999 had statistically significantly 
higher serum 25{OH)D concentrations compared to carriers of all 
risk alleles. 
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